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Abstract 

The Kondo divergences owing to interaction of current carriers with local moments in highly correlated 
electron systems are considered within the Hubbard and s-d exchange models with infinitely strong on- 
site interaction, the many-electron Hubbard representation being used. The picture of density of states 
containing a peak at the Fermi level is obtained. Various forms of the self-consistent approximation are 
used. The problem of the violation of analytical properties of the Green's function is discussed. Smearing 
of the "Kondo" peak owing to spin dynamics and finite temperatures is investigated. 

PACS: 71.10.Fd Lattice fermion models (Hubbard model, etc.), 71. 27. -fa Strongly correlated electron 
systems; heavy fermions, 71.28.-|-d Narrow-band systems; intermediate-valence solids 

The problem of strong correlations and magnetism in many-electron systems is one of the most important 
in the soUd state theory. Since the Hubbard's works of 60's, a great progress has been achieved in under- 
standing electronic structure of systems with strong on-site interaction. Last time, the role of the Kondo 
effect has been discussed within the large-d approach {d is space dimensionality) which reduces the initial 
periodic Hubbard model to an effective Anderson impurity model |Q, ^. Besides the Hubbard bands, an 
important role in the formation of density of states (DOS) picture belongs to a peak at the Fermi level, which 
was found for both half-filled and doped case (the latter case is considered in ^ ) . It should be noted that 
this approach meets with a number of computational difficulties (e.g., consideration of finite temperatures is 
needed, and the low-temperature limit is non-trivial). The structure of the spectrmn in large-d approaches 
is confirmed by the Monte-Carlo calculations. On the other hand, this feature was not reproduced by most 
preceding analytical approaches. In particular, the Hubbard-III approximation Q does not take into account 
contributions of Fermi-like excitations in a proper way because of its single-site character. A detailed anal- 
ysis of this approximation was performed in Refs. ^ within the large-z (2 is nearest-neighbor number) 
expansion. 

In the present paper we present a treatment that is based on the method of equations of motion for 
the many-electron Hubbard operators |0, H and is much more simple than the large-d approach. As the 
zero order this approach reduces to the simplest Hubbard-I approximation. General expressions for 1/z- 
corrections were obtained in Ref. 1^. Unfortunately, the terms with the one-particle occupation numbers 
(which just describe the Kondo effect) were neglected in Refs. ^, and only a classical approximation was 
considered in ReL ^. 

We start from the s-d exchange model with the large s-d coupling parameter |/|, 

W = ^ tkcl^Cka - SiCTcra'cl^Cia' + Ti-d, (1) 

where is the band energy, TCd is the Heisenberg Hamiltonian of the localized-spin system, <t are the Pauli 
matrices. In the limit |/| — > aoo (here and hereafter a = sign/ — ±), it is convenient to pass to the 
atomic representation of the Hubbard operators — where Hsd (the second term in (|^)) takes 

the diagonal form For the electron concentration n < 1, after performing the procedure of projection 
onto the corresponding state space, the one-electron Fermi operators c|g. are replaced by the many-electron 
operators gl^^i ■ These are expressed in terms of the X-operators as 



^ - ^{{S + <jM + 1)/{2S +l)y/^X,{M + ^,+;M), 



M 
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5.1- = ^a{(5-aM)/(25+l)}i/2x,(Af + |,-;M), 

M 

where \M) are the empty states and \ma) are the singly-occupied states with the total on-site spin S + a/2 
and its projection m (which survive at |/| — » aoo), cr = ±. The Hamiltonian of the s-d exchange interaction 
yields a constant energy shift only and can be omitted, so that we obtain 

W = X! ^ksLa^kCTQ + Wd, a = sign/. (2) 

For n > 1 we have to pass to the "hole" representation by introducing new localized spins S = S ± 1/2, and 
the Hamiltonian takes the same form (||) with the replacement tk ^ ~t\i{2S + 1)/(2S' + 1). For S = 1/2, 
I ^ —oo the Hamiltonian coincides with that of the Hubbard model in the case of the infinite on-site 
repulsion U oo with the replacement ik/2 and the almost half-filled band, so that n ^ 5 with S the 
hole concentration. Thus we do not need to discuss the Hubbard model separately (physically, both models 
describe the local- moment situation). 

We calculate the one-particle Green's function for ferromagnetic (FM) state in the energy representation, 

GkaciE) = ({gkaa\gl^J)E- (3) 

Using the commutation relations for for the operators gia-a and spin operators (cf. Refs. ^) we obtain 
the equation of motion 

(£^-ik-q.a)Gk.a(£;) = P.„ + ^^^^ik-q{((<^[^5totq- ^5]^q(^«;^«)]5k-q.a|5La))£ 

q AI 

+ {{Stotq9k-q,-aa\9Lj)E}, (4) 

with S^^, S'j^Qj being the total spin operators (including the contributions of empty and singly-occupied 
states, see for details [^), 

_ 5+l/2 + ac7(g-,t)-W2 

is the total on-site average magnetization. The "Kondo" term comes from the last Green's function in 
the right-hand side of (^, which describes spin-flip processes. Performing decoupling of the next equation 
of motion to first order in the nearest- neighbor number l/z we derive 

{E - tk-q<Ta)((5'to^qffk-q,-<Ta|5La))B = 25 + 1 ^^"^^ ^ ^'^^ ^ l)"k-q-aa](l + tkGk<Ta(£^)) 

-ank-q-aatk-qGwaaiE). (5) 

Here we have neglected TCd, 

Xqo- = ('5'tot-q'5'totq)j '^kiTQ = (.^ko-affkcra ) ■ 

To lowest order (i.e. in the Hubbard-I approximation) we have 

^kcTQ -^cro y (^kCTO ) 

where f{E) is the Fermi function. After substituting (^) and a similar equation for the longitudinal contri- 
bution in (^) (which contains Xq — ('S'tot-q'S'totq)) into (^) we obtain 

bkaa{E) {E)tu' 
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where 



flkcra {E) = Paa + ^ 



q 

t 



+ {2S + l)nk-q- 



q 



(25+1)2 

2 

k-q n 

2S + 1 E — ik-q-tra 



E ^k — q— CTa ^k— q(TQ 



(7) 



The longitudinal fluctuations can be neglected in the saturated FM region |T^, but are important in the 
paramagnetic region (the corresponding results are obtained from the above formulas at (S^^^) = 0) which 
is considered hereafter. The equation for the chemical potential reads 



To smear the logarithmic DOS singularity at the Fermi level, one has to take into account spin dynamics 
which is determined by the Haniiltonian TCd (cf. [pT[ ). To this end, we introduce the normalized spectral 
spin function 

where N-Qii^) is the Bose function. In the far paramagnetic region (where spin correlations can be neglected) 
we can use the simplest spin-diffusion approximation with 



, , 1 Dg2 



where D is the spin diffusion constant (in fact, the results depend weakly on the concrete form of spin 
dynamics). Then we have to put in (|^) 

,J^^ U , ST' f ^ f \ *k-q X + (25" + l)nk-qa 

bi,c.{E) = E-Y fdivKJu;)^^ "^^^^ . (11) 

To simplify numerical calculations, we average the spectral function (^) in q, 

ifq(c.)^ir(c.)=^Xq(^), 

q 

which is sufficient to obtain qualitatively valid results. Indeed, this approximation (which is in spirit of the 
large-d or large-z expansion) reproduces correctly the low-frequency behavior of spin fluctuations which is 
important near the Fermi level. Then a{E) and h{E) do not depend on k. We have taken below D = 0.7c|t| 
{t is the transfer integral). 

In the far paramagnetic region we have 

Then we obtain for 5 = 1 /2 

_ ( 3(1 - n)/4 + 2n, a - 
^^1 3(l-n)/4, a: 
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The results of numerical calculation of the single-particle density of states 

Na,{E)^--lmy^Gk^{E) (12) 

TT ^ — ' 
k 

for the semielliptic bare DOS are shown in Fig. |^. One can see that pronounced density of states peaks 
occur at the Fermi level. The peaks are smeared by both including spin dynamics and finite temperatures. 
It should be noted that, unlike the large-d approach the limit of zero temperature makes no difficulties. 

As well as for small-/ Kondo problem, the "Kondo effect" is connected with the Fermi functions (which 
were not treated in Ref. j^) and is a quantum effect that is small in 1/5. Unlike the FM case, where the 
singularity has one-sided form and non-quasiparticle (incoherent) contributions to the density of state plays 
main role |l^, the logarithmic contribution is symmetric with respect to the Fermi level. Note that after 
formal expansion in I/2: the contributions to the Green's function with the Fermi functions are canceled in 



the paramagnetic phase to first order |12 . 

As discussed in Ref. [|], the approximation (|^). (||) leads to some formal difficulties connected with 
occurrence of an additional singularity of the Green's function in the complex plane. This can result in the 
violation of analytical properties, in particular, of the normalization condition 

L^{{g^aa,9lJ)^Pa (13) 

where 

+00 



i= y" dENa{E). 



The analytical properties turn out to be different for the cases of positive and negative /. For / > the 
singularity lies in the upper half-plane, and for / < in the lower one, so that the normalization condition is 
violated for / < 0, although this violation is numerically not too large (see Table This difficulty does not 
seem to lead to unphysical conclusions in our simple approximation, since the results for the DOS picture in 
both the cases are qualitatively similar. However, it should be stressed that this problem has a rather general 
character and is typical for most calculations which use the many-electron Hubbard representation (or related 
slave-boson and slave- fermion representations which include constraint conditions). It is connected with the 
k-dependence of "perturbation" (band energy it), which changes its sign in the band, and a complicated 
structure of the Green's fucntion (the simple Dyson equation is not valid). By these reasons, the signs of the 
corresponding imaginary parts are not fixed; the problems increase when using self-consistent procedures. 
In particular, such a difficulty should occur in the non-crossing approximation (NCA) for the Anderson 
model where an expansion in the hybridization Vk is constructed. Unfortunately, this problem is usually 
misregarded since the normalization condition is practically never verified. 

The approximation (^), (0), (||) has an unpleasant physical drawback: DOS has the "Van Hove" singu- 
larities at the energies, corresponding to the edges of the Hubbard-I band [^. To remove this drawback, we 
renormalize self-consistently the bandwidths in the resolvents (i.e., in the denominators in (|l0|), (|lT|)) and in 
the Fermi functions by replacing 

E - Pat).-q -^E- P^t^^q. (14) 

Due to this self-consistency procedure, the band edges in the resolvents coincide with those for the total 
Green's function (^). The corresponding numerical results are shown in Fig. ^The analytical properties 
for / > are not violated (Table |^). This approximation is in spirit of large- iV expansion and retains the 
quasiparticle picture, unlike the Hubbard-III aproximation and the self-consistent approximation considered 
below which yield a strongly incoherent behavior. 

Finally, we discuss also the "true" self-consistent approximation where the exact Green's functions are 



included into resolvents. Then we have to replace in the denominators of ( |10[) , (|ll| ) 

E — -Pa^k-q bl^-qa{E) — Ok-qa (-E)tk-q (15) 

and to make the corresponding renormalization of the distribution functions 



rika ^-^J dEf{E)ImGi,c.iE). 
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One can see that the DOS picture (Fig. ^ becomes somewhat smeared even in the absence of spin dynamics, 
but the smearing is not so strong as in the FM case . In this self-consistent approximation, the analytical 
properties are violated for both / > and / < (Table |^). 

To conclude, we have analyzed formation of the "Kondo" peak in narrow-band systems at the Fermi level 
within a simple analytical approach. The numerical results demonstrate sharp energy dependence of the 
density of states near Ep. As well as in the standard Kondo problem, one can expect strong temperature 
dependences of thermodynamic and transport properties (e.g., an enhancement of the effective electron 
mass). This problem needs further investigations. In particular, summation of higher-order perturbation 
corrections with the use of slave boson and fermion representations would be of interest. 

We are grateful to A.O.Anokhin for useful discussions. The research described was supported in part by 
the Grant 99-02-16279 from the Russian Basic Research Foundation. 
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Figure and table captions 



Table 1: Values of the normalization factors L/Pa (see Eq. (|l3| )) for semielliptic bare DOS in (I) the 

approximation (|8|); (II) the self-consistent approximation ([14|); (III) the self-consistent approximation 

(II)- 

an L/ Pa 
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.000 
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.002 
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.000 
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.000 
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.001 


- 0, 


.00 
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.510 


1, 


.395 
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.325 


- 0, 


.02 
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.489 
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.328 


- 0, 


.05 
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.448 
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.337 
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.277 
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.15 
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- 0, 
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.289 


1, 


.192 
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.177 
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Figure 1: The density of states picture for the Green's function and semielliptic bare DOS with S — 1/2, 
n = 0.15 in the absence of spin dynamics for different temperatures (Eqs. (|^), (||), hnes 1, 2, 3 correspond 
to T — 0, 0.01, 0.03) and with account of spin dynamics at zero temperature (Eqs. ( 
a = + and (b) a = — . Inset shows the region near the Fermi energy which corresponds to the peak top. 
Energy and temperature is measured in units of bare half-bandwidth. 



10), (0), line 4); (a) 
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